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Abstract 

A new approach proposed recently by author for the calculation of Green functions in quantum 
field theory and quantum mechanics is briefly reviewed. The method is applied to nonperturbative 
calculations for anharmonic oscillator, (^^-theory, quantum electrodynamics and other models. 



1 Toy model 

To illustrate the general properties of the calculation scheme, which was proposed recently in works 



0^ 



we consider a toy problem of the calculation of " n-particle Green functions" in zero-dimensional theory 
with 0^-interaction, i.e., a problem of calculating the following quantities 



^ ■ G„ = gn/go, (1) 

' where 

/■°° m'^ 

in; g„= / d0</.2"exp{~— 02-A0n- (2) 

J -oo ^ 

^ , A generating function for these quantities is 

cn 

g : gix) - J d0exp{-^02 „ ^ (3) 



At that gn = U=o- 

The generating function g(x) satisfies the differential equation 

, , 2 r, ^d.g 



4A-^ + K-2.)-^-,^0, (4) 

1^ i which is the Schwinger-Dyson equation (SDK) for this toy model. 

A solution of this equation can be written as a coupling constant perturbation (CCP) series g{x) = 
9pert + • ■ ■ + 5peri + ' ' ' whcrc leading ("free field" ) approximation is g'^J^^ — exp{ — ^ log — 2a;|} (up to 

' a multiplier which inessential for calculating of G^*-*) and gp^^t ^ solution of the perturbation iteration 

I scheme equation 

/ 2 n \ ^ypert (i) . , " ilpert 

(m -2x)— — -g^,;, ^-4A / . 



From the point of view of differential equation theory the CCP is attributed to the type of so-called 
singular perturbations. This fact defines, to a considerable extent, both the poor convergence properties 
(the CCP series is an asymptotic expansion at best) and the limited nature of the field of its applicability. 
(In the model under consideration the CCP theory gives good results in the region A < 0.01m'' and nothing 
more.) A perturbation is named to be singular if it contains a higher derivative term. The perturbation 
theory over A is singular in the above sense, since the leading approximation consists in the neglecting of 
the higher derivative term. (See |^ for more discussion.) 

An alternative for the perturbation theory over A can be other iterative scheme that based on an 
approximation of eq. (1) near the point a; = by an equation with constant coefficients. Take as a 
leading approximation the equation 



V — r-^ 1- TO — ; ' = 0. 
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The term Sx^^ will be considered as a perturbation. The iteration scheme will consist in step- to step 
solutions of inhomogeneous equations 

Leading approximation equation has the solution g^^^ = e"^, where a's are roots of the characteristic 
equation 

4Aq!^ + m^a - 1. 

From a condition of the convergence of integral (3) at A ^ we choose the root 

a = ^ = l/m^ + OiX). 

A solution of nth-step equation is looked for in the form (7^*^ = p^^^g^'^^ Taking into account the leading 
approximation equation (5), we obtain the equation for p^'^ 

It is evident from this equation that p^'^ is a polynomial of degree 2i in a;. 

The question about a small parameter for the expansion defined by cqs. (5)-(7) arises. There is no 
manifest small parameter for this expansion , but it is clear, that the expansion approximates well the 
exact solution not only for small values of A. Really, in the strong coupling region A — > 00 the first- 
step one-particle function Gi ("propagator") approximates the exact result with accuracy of 26%, and 
second-step one with accuracy 7%. 

To be more exact, the question about the small parameter should be replaced by the question about a 
convergence of the expansion. But the convergence of this iteration scheme can be easily proved. Notice, 
the iteration scheme defined by eqs. (5)- (7) is equivalent to iterations of the second kind Volterra equation 
with the continuous kernel. The convergence of the iteration of this equation is fulfilled by the textbook 
theorem. 

So, this regular expansion possesses the good convergency properties in the framework of this simple 
zero-dimensional model and, more importantly, is a nonperturbative method of calculations of Green 
functions. 



2 0^-theory 



Let us go to the field theory. Consider the theory of a scalar field (j){x) in the Euclidean space E4 with 
the action 

A{cf>) = I dx{^{d^cf>r + ^ ^^4} (8) 

and with the generating functional of Green functions (vacuum expectation values) 

G{r]) = j D<l>exp{-A + ct)r](l>}. (9) 

Here •q{x^y) is a bilocal source, (prjcf) = j dxdy(j){x)ri(x,y)(j){y). A normalization constant is omitted. The 
nth derivative of G over r] with the source being switched off is the 2n-point (n-particle) Green function. 
The Schwinger- Dyson equation (SDE) for the generating functional Girj) is a corollary of the identity 

Q = j D^j^{ct>{y)eM-A + ctyn<l>})- (10) 



Taking into account the above definitions we get the SDE for (/)^-theory 
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At d = (zero-dimensional theory, or " single-mode approximation" ) the functional derivatives transform 
into usual ones, and eq. (11), after obvious redesignations, reduces to the ordinary differential eq. (4). 
At c? = 1 the model corresponds to the quantum-mechanical anharmonic oscillator. At d > 2 (field 
theory) for the cancellation of ultraviolet divergences the appropriate counterterms should be included 
in the action. The SDE for the theory with counterterms has the form of eq. (11) with the substitution 
A X + 5X, m? rrP' + Sm?, (l + (5z)9^, where 5\,5m? and 5z are correspondingly counterterms 

of coupling, mass and wave function renormalizations. 

Let apply to SDE (11) the same idea about the approximation by an equation with "constant" (i.e., 
independent from 77) coefficients. As the leading approximation equation we will consider the equation 

4A^!^ + (^2 _ 92)^ _ g(o) ^ 0^ (12) 
orfQ'q drj 

and the term 277^ (that contains the source t] explicitly) should be treated as a perturbation. Since 
Green functions are the derivatives of G(ry) in zero and the only behaviour of G near ry = is essential 
, such an approximation seems to be acceptable. The iteration procedure for the generating functional 
G = G'^^^ + G'^"'^' -|- • • • -I- G^*'' consists in the step-to-step solution of the equations 

,/^^^^^.a^^'J:^^a"' = 2,'-^. (13) 

077077 077 07/ 

The solution of the leading approximation equation (10) is the functional 

G(°) = exp{ / dxdyr]{y, x) A^") {x - y)}, (14) 



where A^"^ is a solution of the "characteristic" equation 

4A A(°) (0) A(°) ix-y) + (m^ - 9 2) A(°) {x ~ y) ^ Six - y). (15) 

At d > 2 the quantity A(°^(0) must be considered as some regularization. 

Equation (15) is similar to the equation for the propagator in the leading approximation of the 1/A'^- 
expansion. Certainly, the similarity is completely superficial, since the principle of the construction of 
the approximation scheme is different. 

The solution of equation (15) is the free propagator A^^^ — {^^ — d^)^^ with the renormalized mass 
fj,^ — 77i^-|-4AA(°) (0). The quantity A(°)(0) is defined from the self-consistency condition. The propagator 
is the first derivative of G(7;) over the source t] : Gi = A = ^ |i7=o- As can be easily seen, it is simply 

A^*^-* for the leading approximation. 

Notice, that the generating functional (14) of leading approximation does not possess the complete 
Bose-symmetry. Really, as follows from the definition of generating functional, the Bose-symmetry implies 
on full generating functional the condition 

(16) 



Sr]{y, x)Sr]{y', x') 5r]{y' , x)Sri{y, x') ' 

Evidently condition ( [T^ ) does not fulfilled for G*-"-* defined by eq.(14). The violation of this condition leads 
particularly to the violation of connected structure of the leading approximation two-particle (four-point) 
Green function. 

Such a situation is rather typical for nonperturbative calculational schemes with bilocal source (for 
example, for 1/iV-expansion in the bilocal source formalism), but discrepancy of such type are not an 
obstacle for using these iteration schemes. Really, condition (^ should be satisfied by the full generating 
functional G which is an exact solution of SDE. It is clear that an approximate solution may do not 
possess all properties of an exact one. In given case we have just the same situation. Properties of 
connectivity and Bose-symmetry of higher Green functions, which are not fulfilled at first steps of the 
iteration scheme, "improves" at subsequent steps. For example, the structure of disconnected part of the 
two-particle function is reconstructed as early as at the first step of the iteration scheme. At subsequent 
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steps the correct connected structure of many-electron functions and other corollaries of Bose-symmetry 
are reconstructed. Such stepwise reconstruction of exact solution properties is very natural for the given 
iteration scheme as it is based on an approximation of the generating functional G(r/) in vicinity of 
zero. The Green functions are coefficients of the generating functional expansion in the vicinity of zero, 
therefore only the lowest functions are well-described at first steps of the approximation - at the leading 
approximation the propagator only. Higher many-particle functions come into the play later, at following 
steps, and relation is fulfilled more and more exactly when we go toward exact solution. 

In the general case, the solution of equation for the z-th step of the iteration scheme is the functional 
= PWG^"), where P« is a polynomial in 77 of a degree 2i. Therefore at the j-th step the computation 
of Green functions reduces to solving a system of 2i linear integral equations. 

A solution of the first step equation is G^^^ = F^^^G^^^ where P^^^ = ^Fr]"^ + A'-^^rj. The function F is 
two-particle (four-point) function of the first step, and A*-^^ is the first-step correction to the propagator. 
Eq.(13) at I = 1 gives us a system of equations for F and A*^^^. Equations for F and A*^^^ are simple 
linear integral equations. The exact form of solutions of these equations see in §. At A ^ the first 
step propagator reproduces correctly the first term of the usual CCP theory. 

At d = 1 the model with action (6) describes the quantum-mechanical anharmonic oscillator. Ultra- 
violet divergences are absent, quantities of A'^*'^(0) type are finite and the above formulae are applied 
directly for the computation of Green functions. 

To calculate a ground state energy E one can use the well-known formula 

dE 

— = G2(0, 0,0,0), 

where G2 is the four-point (or two-particle) function. Integrating the formula with a boundary condition 
^^(A = 0) = m/2 taken into account, one can calculate the ground state energy for all values of the 
coupling (see [||). 

At A — > the first step calculation reproduces the perturbation theory up to the second order. 
At A ^ 00 : E = eoA^/^ + 0{X~^^^), and eq = 0.756. The coefficient eo differs by 13% from the exact 
numerical one eg^'^'=* = 0.668. At A/m'^ — 0.1 the result of the calculation differs from the exact numerical 
one by 0.8% and at X/m^ — 1 differs by 6.3%. Therefore, the first step calculations approximate the 
ground state energy for all values of A with the accuracy that varies smoothly from (at A 0) to 13% 
(at A 00). 

At d > 2 action (6) should be added by counterterms for the elimination of ultraviolet divergences. 
There is no need to add a counterterm of wave function renormalization for the leading approximation, 
and the equation of the leading approximation will be 

x2^(0) Ar'(o) 
4(A + SXo)-^ + (Sml + m^- 9^)^ - G^ = 0. (17) 
077077 077 

At i > 1 the counterterms SXi, Smf and Szi should be considered as perturbations. Therefore, the 
corresponding terms should be added to the r.h.s. of equation (11). So, the first step equation will be 

4(A + ao)V^ + ('5™o + "i'-a')^-GW = 
077077 077 

= 277— Sml —c — + ^zid^ —c 4(5Ai . (18) 

077 077 077 077077 

For the super-renormalizable theory (d = 2 and d = 3) it is sufficient to add counterterms of mass 
renormalization and wave function renormalization, i.e. SXt = for all n. The normalization condition 
on the physical renormalized mass /i^ gives us a counterterm of the mass renormalization in the leading 
approximation. This counterterm diverges logarithmically at d = 2 and linearly at d = 3. The countert- 
erm Szi is finite at d = 2,3. The counterterm Smf diverges as that of the leading approximation does, 
namely, logarithmically at d = 2 and linearly at d = 3. 

At d = 4 besides the renormalizations of the mass and the wave function a coupling renormalization 
is necessary. Due to the presence of the counterterm SX the normalization condition on the renormalized 
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mass /i^ for the leading approximation becomes the connection between counterterms Sm^ and SXq. 
Counterterm i5Ao (and, consequently, Sttiq) will be fixed at the following step of the iteration scheme. 

A solution of the equation for the four-point function at d = 4 diverges logarithmically, and a 
renormalization of the coupling is necessary. The equation for F contains the counterterm SXq only. 
Therefore by defining a renormalized coupling Ar as a value of the amplitude in a normalization point 
we obtain the counterterm of the coupling renormalization SXq and the renormalized amplitude. Taking 
the renormalization of the two-particle amplitude in such a manner, one can solve the equation for A*^^^ 
and renormalize the mass operator in correspondence with the general principle of normalization on 
the physical mass. But in four-dimensional case one gets an essential obstacle. At the regularization 
removing, ^Aq — > —A, and the coefficient A -I- SXq in the leading approximation equation (17) vanishes. 
The same is true for all the subsequent iterations. The theory is trivialized. One can object that an 
expression 

is really an indefinite quantity of • oo type, and the renormalization is, in the essence, a definition of 
the quantity. But it does not save a situation in this case since the renormalized amplitude possesses 
a nonphysical singularity in a deep-euclidean region (it is a well-known Landau pole). The unique 
noncontradictory possibility is a choice Ar ^ at the regularization removed. This is the triviality of the 
theory again. The triviality appears almost inevitably in an investigation of (/)|-theory beyond the CCP 
theory and is a practically rigorous result. Notice, that contrary to the CCP theory which is absolutely 
nonsensitive to the triviality of the theory, the method proposed leads to the triviality already at the first 
step. 



3 Quantum Electrodynamics 

The Lagrangian of Quantum Electrodynamics (QED) in Minkowski space-time with a gauge fixing term 
has the form 

C - -^F,,,Ff,, - -L(a^A^)2 + ^{id - m + eA)^. (20) 
4 Zdi 

Here F^^ — d^A^, — di,A^, A = A^7^, i/i = "0*70, 'ti is an electron mass, e is a charge (coupling 
constant), di is a gauge parameter, 7^ are Dirac matrices. For notation simplicity we write all vector 
indices as low ones. 

A generating functional of Green functions is 

G(J,7?) = j i?(^,^,A)exp*{ j dx(c + J^{x)A^{x)) - J dxdy^pf'iy)7j^^{y,x)i:"ix)}. (21) 

Here J/i(a:) is a source of the gauge field, and ri^°'{y,x) is a bilocal source of the spinor field (a and 
f3 are spinor indices). Normalization constant omitted. 

Functional derivatives of G with respect to sources are vacuum expectation values 

SG . „ „ SG 



SDKs for the generating functional of Green functions of QED are 

S{x - y)G + {rd m)-^ + -7, ., /'f, r " / dx'rjix, x')-^^ = 0. (24) 

dr-i{y,x) I dJf_,{x)dri{y,x) J dr]{y,x') 

(Here and everywhere below 9^ denote a differentiation with respect to variable x.) Let us resolve SDE 
(p3|) with regard to the first derivative of the generating functional with respect J^: 



1 SG 



i 5Jf^{x) 



= - [ dxiD" (x - Xi){j.{xi)G + ietr-f.—^ — -}, 



(25) 
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and substitute it into the second SDE ([2^). As a result we obtain the "integrated over A^" (in the 
functional-integral terminology) equation 

d(x - y)G + [id - m)— + — / dxiD [x - xi)-/,,— r tr7^— = 

5ij[y,x) I J ^ '^dT^[y,x) d'q{xi,xi) 

{SG SG 
Vix.xi)— + eD" {x - xi)J^{xi)jf,— -I. (26) 
dT][y,xi) dr]{y,x)) 

Exploiting Fermi-symmetry condition 

S^G 6^G 



Sr]l^"'{y, x)5r]P'°'' (y', x') S-qf^'"" {y', x)6f]'^°'' {y, x') ' 
let us transform cq.(p9) in following manner: 



(27) 



S{x - y)G + {id - to) + ie^ [ dxiD" {x ~ xi)^^—^ -^y- 



Sri{y,x) J '"^ Sr]{xi,x) Sri{y,x 



= [ dxi{T]{x,xi) + eD" {x - xi)J^{xi)-f^ \. (28) 

From the point of view of exact solutions equations (^ ) and (^) are fully equivalent since the transition 
from eq.(p6|) to eq. (p8|) is, in essence, an identical transformation. However, it is not the case for the 
used iteration scheme since, as for Bose-symmetry condition (|l6| ) as, Fermi-symmetry condition (2^) is 



fulfilled only approximately at any finite step of the iteration scheme. Therefore, eqs. (|6|) and (^) lead 
to different expansions. Eq.(p6|) gives to the calculational scheme, which on the language of Feynman 
diagrams of perturbation theory is analog of the summation of chain diagrams with fermion loop. This 
version is named "calculations over perturbative vacuum" since a unique connected Green function of 
the leading (vacuum) approximation is the free electron propagator. This scheme leads, as for 0^-theory 
as (see above, section 2), to Landau pole and triviality at the first step of the iteration scheme (see Q 
for more details). 

The second version of the iteration scheme, which is based on eq. (|2^) , gives us a fruitful and " insen- 
sitive to triviality" scheme of calculation of physical quantities. This version is named "calculations over 
nonperturbative vacuum" since the electron propagator of the leading vacuum approximation is a solution 
of a non-trivial nonlinear equation. For this scheme a calculation of two first terms of expansion of the 
vertex function in photon momentum for chiral-symmetric vacuum have been performed in work . This 
calculation has allowed to obtain a simple formula for anomalous magnetic moment: /2 = a/(27r — a), 
where a is the fine structure constant. Also, for a linearized version of the theory (see, for instance, Q) 
the problem of dynamical chiral symmetry breaking have been investigated in work . The calculations 
are performed for renormalized theory in Minkowski space. In the strong coupling region a > 7r/3 the 
results correspond to earlier investigations performed in Euclidean theory with cutoff (see |6j): solutions 
arise with breakdown of chiral symmetry. But for the renormalized theory a solution with breakdown of 
chiral symmetry is also possible in the weak coupling region a < tt/3 with a subsidiary condition on the 
value of a which follows from the gauge invariance (see Q for more details). 



4 Other models 

Some other models have been investigated by proposed method in works and (7|-|P. Famous 

Gross-Neveu model with the Lagrangian 

£ = ^,*aV, + ^(V5,V',)' (29) 

have been investigated by this method in work at D = 2, 3, 4 and finite N, where D is space dimension 
and A'^ is a number of flavors. The results were following: a spontaneous symmetry breaking is shown 
to exist in £) = 2, 3 and the running coupling constant is calculated. The four dimensional theory turns 
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seems be trivial. These results exhibit the efhciency of the method and are the finite N generalization of 
the known results obtained in the framework of expansion. 

The greatest interest from the physical point of view presents applications of the method to study 
gauge theories in nonperturbative region. First steps in this direction were made in works In 
work a generalization of the Higgs mechanism which takes into account the contributions of gauge 
field vacuum configuration into the formation of the physical vacuum was considered. For the Abelian 
Higgs model the triviality bound niH < 1.15m^ was found. In works a non-Abelian SU{2) gauge 

theory was considered, and a mechanism for the dynamical mass generation of a non-Abelian gauge 
field which was based on taking into account the contributions of the gauge field vacuum configurations 
into the formation of the physical vacuum was proposed. These investigations are needed in a following 
elaboration. A most winning field of application of proposed method seems to be a problem of dynamical 
chiral symmetry breaking in gauge theories. 
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